Properties of the general NHDM. I. The orbit space 
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Abstract 

We study the scalar sector of the general A^-Higgs-doublet model via geometric con- 
structions in the space of gauge orbits. We give a detailed description of the shape of 
the orbit space both for general N and, in more detail, for = 3. We also comment on 
Oh| remarkable analogies between NHDM and quantum information theory. 

1 Introduction 

> , 

^ . The Standard Model relies on the Higgs mechanism to realize the electroweak symmetry 
! breaking (EWSB). Many variants of the Higgs mechanism have already been suggested, but 
it is not presently known what particular variant is realized in Nature. 
^ ■ While waiting for the LHC to give the first hints of the Nature's choice, theorists must 

be duly prepared to safely interpret the future LHC data. This implies, in particular, that 
theorists must be aware of all essential possibilities which can be realized within a chosen 
model for EWSB. Since non-minimal Higgs sector usually involve many free parameters, it 
^ , is highly desirable to analyze the chosen model in its most generic formulation, allowing 
^ ' for all possible degrees of freedom. This general analysis of a specific model should show 
^ ■ which phenomenological consequences are universal and which are sensitive to values of the 
parameters, which symmetries can in principle arise in the model and how they are broken, 
which properties hold only at the tree-level and which survive the perturbation series. When 
this general structure of a model is well understood, one should proceed further and restrict 
the model by taking into account existing experimental constraints. 

Unfortunately, such an exhaustive analysis is hardly feasible for many non-minimal Higgs 
sectors. A very representative case is given by the 2HDM, [Il|2l|3]. Here, the straightforward 
algebra fails even at the very first step, because the Higgs potential cannot be minimized 
explicitly in the general case. As a result, for a long time only relatively simple variants of 
the 2HDM were analyzed, while the most general 2HDM remained barely studied. In the last 
several years a number of tools were developed which led to many insights into the properties of 
the general 2HDM. These methods were based on the idea of the reparametrization symmetry, 
or basis-invariance, of the model: a unitary transformation between the Higgs doublets changes 
the parameters of the model, but nevertheless leads to the same physical properties of the 
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observable particles. This idea can be implemented via the tensorial formalism at the level of 
Higgs fields [H O El |7] or via geometric constructions in the space of gauge-invariant bilinears 
[El [9l [ini im US]- In the latter case the formalism was extended to include non- unitary 
reparametrization transformations [13l|14l[l5], which revealed interesting geometric properties 
of the 2HDM in the orbit space equipped with the Minkowskian metric. 

It is a natural idea to extend these successful techniques to N doublets. The general 
NHDM is obviously more involved than 2HDM, both at the level of scalar sector and Yukawa 
interactions (see examples in [T6l[T71[l8]). Some properties of the general NHDM were analyzed 
in [iniEOllIOlEHEHlEl], with a special emphasis on CP- violation, [241 E]. However, a method 
to systematically explore all the possibilities offered with doublets was still missing. 

In principle, generalization from 2HDM to NHDM is straightforward in the tensorial for- 
malism, however it is very difficult to translate tensorial invariants into physical observables. 
On the other hand, the geometric approach in the space of bilinears offers a more appealing 
treatment of the Higgs potential, but the shape of the NHDM orbit space is rather compli- 
cated and has not been fully characterized so far. In this paper we fill this gap by studying in 
detail the algebraic and geometric properties of the NHDM orbit space. Many of these results 
will be used in the companion paper [25j where the minimization problem and the symmetry 
breaking patterns of the Higgs potential of the general NHDM are analyzed. 

The paper is organized as follows. Section 2 is devoted to three distinct, but interrelated 
approaches to description of gauge orbits in the space of Higgs fields. Then, in Section 3 
we construct the orbit space of NHDM as a certain algebraic manifold and discuss at length 
its algebraic and geometric properties. In Section 4 we treat the specific case of 3HDM in 
even greater detail, aiming not only at a concise algebraic description of the orbit space, but 
also trying to gain an intuitive understanding of its shape. In the final Section we draw our 
conclusions. 

2 Describing Higgs fields in NHDM 

The scalar potential of the NHDM is constructed from gauge-invariant bilinear combinations 
of the Higgs fields. The space of these combinations can be described in three algebraically 
different but closely related ways: via representative Higgs doublets, via a i^-matrix, and via 
a vector in the adjoint space. In this Section we describe and compare these three ways. 

2.1 Field space and gauge orbits 

The scalar content of the general NHDM consists of complex Higgs doublets with elec- 
troweak isospin Y = 1/2: 



The total dimensionality of the space of scalar fields is AN. Since the Higgs lagrangian is EW- 
symmetric, we can perform any simultaneous intradoublet SU{2) x U{1) transformation inside 
all doublets without changing the lagrangian. If we take a generic point in the Higgs space 
and apply all possible EW transformations, we will get a four- dimensional manifold called the 
(gauge) orbit. Thus, the entire 4iV-dimensional space of Higgs fields is naturally "sliced" into 
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non-intersecting orbits. The resulting set of orbits is a (4A^ — 4)-dimensional manifold called 
the orbit space. 

In principle, (pa are operators. However, when minimizing the Higgs potential, we will look 
for vacuum expectation values of the Higgs fields (0a), which are c-numbers. Then, we can 
characterize each orbit by a specific representative point in the Higgs space: 

^-(°)- 

This point (and therefore, the entire orbit) is characterized by AN — 4 real parameters: 
values of Va, N — 1 values of u^, a > 1, — 1 phases ^a, a > 1, and N — 2 phases rja, a > 2. 

It is well known that if at least one Ua 7^ 0, such a point corresponds to the charge- 
breaking vacuum, in which the electroweak symmetry is broken completely and the photon 
acquires mass. If we insist that the vacuum be neutral, we must set all Ua = 0. Thus, the 
representative point of a generic neutral orbit is 

It is characterized by parameters Va and N — 1 phases ^a, making the dimensionality of the 
neutral orbit space equal to 2N — 1, which is 2N — 3 units less than the dimensionality of the 
entire orbit space. 



2.2 Reparametrization freedom 

So far we have mentioned only the electroweak SU{2) x f/(l) transformations between the 
components of each doublet. Since the potential is an EW scalar, such transformation do not 
affect the parameters of the potential. 

Now consider the SU{N)h group of transformations that mix the doublets without affect- 
ing their intradoublet structure (index H stands for "horizontal"). This transformation sends 
a given Higgs potential to another viable Higgs potential with different coefficients. Such a 
transformation is called a reparametrization transformation, or a horizontal space transfor- 
mation, or a Higgs-basis change. The key property of this transformation is that although it 
reparametrizes the potential, it leaves the physical observables invariant [U [5] . This property 
is known as the reparametrization invariance, or the Higgs-basis invariance, of the model. The 
same is true for anti-unitary transformations as well, so one can state that the reparametriza- 
tion group of the general NHDM consists of all unitary and anti-unitary transformations acting 
in the space C^: (pa ~^ Uab4>b and 0^ — )■ Uab4>l- The antiunitary transformations are also known 
as generalized CP transformations [lOl [39] . 

Reparametrization transformations link different gauge orbits. If we pick up a specific point 
in the gauge orbit space, then by applying all possible reparametrization transformations we 
can reach many other points in the orbit space. Thus, the orbit space itself becomes split into 
non-intersecting SU{N)-oThits. We will study this stratification in more detail in Section 13.41 

For any neutral orbit parametrized by Va, C,a according to ([3]), we can find a reparametriza- 
tion transformation that brings it to a "canonical form" (also known as Higgs basis) 
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Obviously, we have equivalent canonical forms depending on which doublet has the non-zero 
value. These equivalent forms can be related to each other by a discrete subgroup (permuta- 
tion) of the reparametrization group. 

Equivalently, any point in the charge-breaking orbit space can be brought to its own 
canonical form 

0i=fM, (j)2=('^], 0a= ( Q j for a>2, t;2 = ^t;2, m2^^w2_ 

a a 

There are iV(A^ — 1)/2 such canonical choices. To avoid double counting of equivalent canonical 
forms within such choices we can restrict v"^ > v?. The neutral orbit space corresponds to 
the limit v? — )■ 0. The other extremum, v? = f^, corresponds to a rather special space of 
^''maximally charge-breaking^^ (MCB) vacua. 



2.3 i^'- matrix formalism 

Following [9l [TTl [10] , we represent the electroweak-scalar bilinears defined above as components 
of a complex hermitean N x N matrix: 

Kab ^ (0l0a) . (6) 

Several important properties of the i^-matrix were proven in [101 El]: 

• It is hermitean positive-semidefinite matrix. 

• Its rank is two for a generic (charge-breaking) vacuum and one for a neutral vacuum. 
This result stems from the fact that we deal with electroweak doublets and not higher 
representations of the gauge group. 

In other words, the i^'-matrix has at most two non-zero (and positive) eigenvalues, while the 
other N — 2 eigenvalues are zero. For the neutral vacuum, one gets only one non-zero (positive) 
eigenvalue and — 1 zeros. The maximally charge-breaking vacua can be defined in terms of 
/T-matrices by [TiK]^ = 2Ti{K^). 

A reparametrization transformation U G SU{N)h acting on doublets (j) transform also the 
fC-matrix according to the adjoint transformation law: 

K UKUl (7) 

In particular, starting from any i^'-matrix, one can always find such a transformation that 
diagonalizes it: 

K = diag(f ^, 0, . . . , 0) for neutral orbit space; 

K = diag(f ^, M^, . . . , 0) for charge-breaking orbit space. (8) 

These diagonal i^- matrices correspond to the "canonical" orbits (jl]) and (E]). 

From the observation that the if-matrix is a hermitean and rank 2 matrix, we can also 
deduce the dimensions of the neutral and charge breaking orbit space. Let us suppose 
rank(i^') = 2. Then among the A^ lines (columns) of K, there are at most two linearly 
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independent. We can arbitrarily choose them to be the first and second hnes (columns) by re- 
labeling the doublets. The remaining lines (columns) can be rewritten as linear combinations 
of lines 1 and 2, and, because of hermiticity, the expansion coefficients are not arbitrary but de- 
termined by the elements of the linearly independent lines. This fact is proved in appendixiBl 
Therefore, we can choose the set 

Kla = a=l,...,iV, 

K2i> = 0l02, & = 2,...,iV, (9) 

as the set of 4N — 4 algebraically independent gauge invariants Kab = 4>\4>a^ noticing that Kn 
and K22 are real while the rest are complex. Thus the dimension of the charge breaking orbit 
space is AN — 4. The neutral orbit space corresponds to the subcase of rank(i^) = 1 where we 
can choose the first line (column) to be the linearly independent line (column), which implies 
that the dimension of the neutral orbit space is 2N — 1. The dimensions of the neutral and 
charge-breaking orbit spaces can be also deduced using the isotropy groups acting on K; see 
details in Section |33] below. 



2.4 Adjoint representation 

Yet another look at the orbit space of the A^-Higgs-doublet model is offered by the adjoint 
representation of the reparametrization group SU{N)h- 

Since the i^T-matrix ([6]) is a hermitean N x N matrix, it can be decomposed as 

K = ro- ^l^^^lM + r,\,, i = l,...,N' -1. (10) 

Here A, are generators of SU{N) satisfying relations 

2 

AiAj = —dijl^ + ifijk^k + dijk^k ■ 

The coefficient in front of the unit matrix in (fTOj) is chosen for future convenience. The values 
of tq and rj can be extracted from the if-matrix: 

^o=V^^'^' = ^Tr[irA.] . (11) 

If the space of i^-matrices were generalized to the space of all N x N hermitean matri- 
ces, i.e., Mh{N), Eq. (fTOl) would define a linear and invertible map between Mh{N) and the 
space of all possible real vectors = (ro,rj), i.e., M.^ [21], which we call the adjoint space. 
The notation r'^ alludes to the Minkowski-space formalism similar to what was developed for 
2HDM in [I3l HH [11], [15] . Since we limit ourselves only to the (anti)unitary reparametrization 
transformations, we are not going to use this formalism here, although we do think that it 
might be useful in NHDM. Here we use r'^ just as a short notation of (ro,rj). 

When Mh{N) is restricted to the space of positive semi definite rank 2 ii"- matrices, which 
we call the K-space, the mapping ( ITOi) from the i^-space to the adjoint space is no longer 
surjective and its image defines a manifold embedded in which will be denoted by V$ [21j 
(orbit space). 
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The shape of V$, which we study in Section |2] below, is rather comphcated. However, 
to cast the first glance at it, let us remind the reader of the situation in 2HDM. There, the 
i^T-matrix was a hermitean 2-by-2 matrix, and the condition ranki^T < 2 was automatically 
satisfied. The positive-semidefiniteness led to 

TrK > , {TrKf - Ti[K^] > , (12) 

which in the adjoint representation translated into 

ro > , rl>7^, (13) 

where r'is just an index-free label for r^. Thus, in 2HDM the orbit space was represented by 
(the surface and the interior of) the forward lightcone in M^. 

With more than two doublets, the rank-2 requirement generates extra conditions to be 
imposed in addition to (fT2|) . These conditions are formulated as equalities (not inequalities) 
among traces of powers of the ii"-matrix up to K'^ , see details in Section [31 In the adjoint 
space they will translate into a set of algebraic equations on components of r'^. The orbit 
space still lies on and inside the forward lightcone ( JT3l) . but it occupies only a certain region 
inside it. 

A unitary reparametrization transformation of (pa keeps tq unchanged but leads to a ro- 
tation of the vector r in An anti-unitary reparametrization transformation adds to 
that a reflection of A^(A^ — l)/2 components of r. Since the map SU{N) — t- SOlN"^ — 1) is 
not surjective, not all rotations of f can be induced by a reparametrization transformation 
of the doublets. Thus, the reparametrization group in the adjoint space is a certain (and 
rather small) subgroup of the full rotation group 0{N'^ — 1). This is also reflected in the 
rather complicated shape of the orbit space itself, which is invariant only under rather special 
rotations. 

As discussed above, the ii'-matrix of any canonical neutral (jl]) or charge-breaking 
orbit is diagonal. Its decomposition fIlOp involves only diagonal matrices Aj, i.e. the Cartan 
subalgebra of su{N) together with the unit matrix. Thus, the canonical orbits correspond in 
adjoint space to a certain (A^ — l)-dimensional section through the root space of su{N), to be 
discussed in Section [3731 

2.5 Higgs potential 

Let us also discuss how the Higgs potential is described within each of these approaches. 
The generic Higgs potential in NHDM can be written in a tensorial form [H [5]: 

V = n,(0i0fe) + z-aU^l^bM^d) , (14) 

where all indices run from 1 to A^. The potential is constructed from A^^ bilinears (0j^0b), 
and therefore it can be viewed as deflned in the orbit spac^. Coefficients in the quadratic 
and quartic parts of the potential are grouped into components of tensors Yab and Zatcd, 
respectively; there are A^^ independent components in Y and N'^[N'^ + l)/2 independent 
components in Z. 

^The idea to switch to the orbit space in order to simphfy the task of a group- invariant potential minimiza- 
tion is rather old, see [55] and references therein. 
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Within the i^-matrix approach, the potential is still based on the same tensors and can be 
written symbolically as 

V = Tt[YK] + TrTr[ZK ® K] , (15) 

where "TrTr" indicates the traces over the two pairs of indices. 

In the adjoint space, one replaces a pair of doublet indices by the index n = (0,z). The 
Higgs potential take the following form: 

V = -M^r'' + ^A^^r^^r" = -{MoTq + M,n) + ^ (Aoor^ + 2Ao^ron + A.,r,r,) . (16) 

The scalar Mq and the vector Mj are essentially Y^b of (HM written in the adjoint space. The 
total number of free parameters in Mq and the vector Mj is 1 + (A^^ — 1) = A^^. The scalar 
Aqo, vector Aqi and symmetric tensor Aij represent Zabsd in the adjoint space; their parameter 
counting gives 1 + (A^^ - 1) + (A^^ - 1)A^V2 = (N^ + l)N^/2. We stress again that we do not 
use Minkowskian metric in this paper, so all the contractions of pairs of indices are understood 
in the Euclidean sense. 

Let us discuss the relation among the three ways of describing the Higgs field configurations 
in NHDM. 

Working in terms of fields parametrized as ([2]) and ([3]), one can easily describe the entire 
orbit space available in terms of Va, Ua and relative phases. The price one pays for this 
simplicity is that the Higgs potential involves tensors Yab and Zabcd^ whose properties are very 
far from being intuitive. For example, even in 2HDM one must resort to long computer-assisted 
algebraic manipulation in order to formulate the explicit CP-invariance of the potential [6]. 
One can expect that understanding the Higgs potential of NHDM will be even harder. 

In the adjoint space the description of the orbit space becomes much more complicated. 
On the other hand, the treatment of the Higgs potential is dramatically simpler. For example, 
in 2HDM one could easily formulate and prove conditions for existence of a symmetry, one 
could derive theorems about the number and coexistence of minima, etc. Remarkably, the 
same treatment can be extended straightforwardly to NHDM, which will be the subject of 
the companion paper [25]. Thus, the orbit space description represents the only essential 
complication on the way to understanding the properties of the scalar sector in generic NHDM. 

The /^-matrix formalism lies somewhere in between. The i^'-space resembles V$ in the 
adjoint space, but the Higgs potential is still written in a non-intuitive tensorial form. However, 
many lines of argumentation can be started at the level of i^-matrix. 

2.6 NHDM as a quantum information-theoretic problem 

It is a remarkable and perhaps an underappreciated fact that the mathematics behind con- 
structing the orbit space in the A^-Higgs-doublet model is very similar to what is studied in 
the Quantum Information Theory (for introduction, see recent textbooks |271|28]) and in the 
so-called geometric quantum mechanics, [29j . 

In quantum information theory one studies quantum evolution of interacting A^-level quan- 
tum states, called qudits (for N = 2 and A^ = 3 one speaks of qubits and qutrits, respectively), 
which are not necessarily isolated from the environment. One of the basic problems here is to 
describe the space of states of a single A^-level qudit. In general, it is described by a hermitean 
positive-semidefinite N x N density matrix p, which satisfies certain axioms, [30| [271 12H] • For 



7 



pure states rank/) = 1, while for a mixed state 1 < rankp < A^. The density matrix can also 
be decomposed in the basis of the algebra su{N), similarly to ffTOl) : 

p=j^lN + PiXi, t = l,...,N^ -I, (17) 

where the coefficient in front of the unit matrix is fixed by condition Trp = 1. The vector p is 
known as the coherence vector, or the Block vector. All possible Bloch vectors occupy a region 
in the (A^^ — l)-dimensional space called the (generalized) Bloch ball. It is remarkable that 
only recently the structure of the Bloch ball was analyzed for N > 2, [31], |32], [331 EH]- 

Many of these objects have counterparts in HNDM. The density matrix corresponds to 
an appropriately normalized i^-matrix. The neutral vacuum of NHDM corresponds to a pure 
states of a qudit, while the charge-breaking vacuum corresponds to a mixed qudit state with 
a rank-2 density matrix. Higher rank density matrices do not have their counterparts in 
NHDM. The coherence vector for the A^-qudit is analogous to the adjoint space vector f, and 
the generalized Bloch sphere is then just another name for the gauge orbit space. 

It is well possible that the analogy between the quantum information theory and NHDM 
can be pursued further. In any case, we believe that the elaborate mathematics of quantum 
theory can generate new insights into the properties of NHDM or similar problems in particle 
physics. 



3 Orbit space of NHDM 

3.1 Geometric description of the orbit space 

Let us start by describing some geometric properties of the orbit space of NHDM V$ in the 
adjoint space of vectors r^, which can be inferred directly from the definitions. 

Let us first note that the orbit space has a conical shape: if point r'^ G V$, then ar^ G V$ 
for any a > 0. Therefore, in order to understand the shape of V<j>, it is sufficient to study its 
tq = const section at any positive tq. To this purpose, we switch to the (A^^ — l)-dimensional 
space of normalized vectors rij = Vi/vo. The neutral orbit space then lies on the surface of the 
unit sphere = 1, while the charge-breaking orbit space occupies a region strictly inside it. 

It is plain to see that a point in the neutral orbit space in the ra-space is parametrized 
by A^ independent complex numbers up to an overall (complex) factor. In other words, the 
points of the neutral orbit space are in one-to-one correspondence with complex rays in 
passing through the origin, which form the (A^ — l)-complex-dimensional complex projective 
space CP^~^. Thus, the neutral orbit space in the n-space has the shape of CP'^"^ embedded 
into R^'-\ 

The entire orbit space of NHDM can be reconstructed from the neutral orbit space by an 
operation, which we call self-join. By definition, the self-join of a set of points S in an affine 
space is a union of points lying on straight line segments drawn between all pairs of points of 
S. Now, let us pick up two points from the neutral orbit space, whose i^-matrices Ki and K2 
are not proportional to each other. Consider the open interval of i^'-matrices lying between 
them: 

K = cKi + {l- c)K2 , c G (0, 1) . (18) 
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Such i^-matrices are necessarily rank-2 matrices. Inversely, any rank-2 i^-matrix can be 
written (not uniquely) as a linear superposition with positive weights of a pair of rank-1 K- 
matrices. Since the map from the i^-space to the orbit space is linear, the same construction 
holds in the r'^-space, which proves that the entire orbit space is a self-join of the neutral orbit 
space. In loose terms, the charge-breaking orbit space is "stretched" on the wireframe of the 
neutral orbit space. 

Note that this construction is similar, but not completely analogous, to the convex hull that 
arises in the quantum information theory, where the density matrices can have an arbitrary 
rank. It also means that the NHDM orbit space does not possess the strict convexity. 

3.2 Algebraic description of the orbit space 

At the level of i^-matrix, the defining criterion is that the i^-matrix is a positive-semidefinite 
matrix with rank < 2 [TOl |2l]. In other words, it requires that there be at most two non-zero 
eigenvalues, which must be positive. Following [311 [32], we write the characteristic equation 
for the i^-matrix as 

N 

det(Al -K) = \^ + ^{-iysj{K)\^-^ . (19) 

The coefficients can be written as products of the eigenvalues Aj, i.e., the roots of the 
characteristic equation, 

n 

l<ii<-<i„<N 3=1 

as well as in terms of traces of powers of the i^'-matrix: 

(_-\\n-l 

s^{K)=^-^—i:i[K^ + y^{-iys,{K)K^-^], n = l,...,iV. (21) 

For example, si{K) = Tr[i^'] and S2{K) = | (Tr^[i^'] — Tr[i^'^]). When written without vari- 
ables, the identification s„ = Sn{K) will be assumed. In general, positive semi-definiteness of 
matrix K is equivalent to non-negative values of all Sn- In our case, the requirement that the 
i^T-matrix has rank < 2 is equivalent to 

Si{K) > , S2{K) > , Sn{K) = for all 2 <n<N. (22) 

Since K is hermitean and hence diagonalizable, the minimal annihilating polynomial is, instead 
of Eq.dH]), 

K[K^-Si{K)K + S2{K)l]=0, (23) 

which automatically guarantees Eq. f l22p . For the neutral orbit space we require that there be 
only one positive eigenvalue, i.e., 

si{K) > , SniK) = , for all 2 < n < , (24) 

which can be summarized by 

K^ = si{K)K. (25) 
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From Eqs. (1191) or (1211) . it is clear that the coefficients s„(K) are functions of K invariant by 
the reparametrization group action in Eq. ([7]). Because of the positive semi-definiteness of K 
and rank(i^') < 2, such action divides the space V<i. C M^' into SU{N) orbits. Each of these 
orbits can be uniquely characterized by the set {si, S2} of SU{N) invariants (reparametrization 
invariants), since the other Sn are null. If K were allowed to be a general hermitian matrix, 
then all the set {sn} of N invariants would be necessary to characterize all the orbits. The one- 
to-one correspondence between an orbit and a set of invariants applies because they uniquely 
define the eigenvalues of the matrix K, in a given order, and the set of all matrices with the 
same eigenvalues are conjugated to the same diagonalized matrix, forming then one orbit. The 
invariants can be calculated and written in terms of in Eq. (ITOl) : 



SiiK) = \/4^ro, (26) 



iV-1 

d2 



S2{K) = rt,-^\ (27) 



2 



s.{K) = |d,^.,r,r,r,--^^ro[r---^], (28) 
s.iK) = -^^^.r,r,r, + + ^yi^^i]'^ - |] , (29) 



where we defined the totally symmetric tensors 



rll...„-^Tr[{A,...A.„K]. (30) 



The symbol { }+ denotes the sum of strings of A's with all possible permutations of indices 

if = and rg) 



ii to in- We can easily identify F^-f = 6ij and vf} = dijk- Notice Eq. (129!) already assumes 



S3 = in Eq. (1281) . Therefore, we can define the orbit space V$ as the set of vectors G 
that satisfy the set of equalities and inequalities of Eq. (12^ , explicitly given by 

r-o>0, r2-r^>0, di^kVirfk + -^j^==^ro{rl - ?,r^) = Q , ■■■ , (31) 

where each equality s„ = 0, n > 3, gives an algebraic equation of order n in r^. A systematic 
procedure to find the the higher order equations s„ = is given in appendix O In general, 
Sn{K) = is equivalent to 

SninX.) - {-irt-^^ - ^0(1 - !)] = , n > 2 , (32) 

where the ffist term depends only on r and it can be written as a sum of terms containing 
contractions of the tensors in Eq. fISU]) up to order n, see appendix O We can see that, for 
each set of values of the invariants tq and r^, Eqs. flHT]) and fl52]) lead to a set of algebraic 
equations of order n < N that defines a manifold on MJ^"^ . Each of these manifolds constitutes 
a single SU{N) orbit in the orbit space because all the available invariants si, S2, are fixed. In 
particular, the neutral orbit space is a particular SU{N) orbit for which the second condition 
of Eq. (13T]) becomes an equality: 

rl-f^ = 0, (33) 



10 



meaning that the neutral orbit space must he on the forward hghtcone. 

For a complete characterization of the orbits, it remains to specify the range of variation 
for {rQ,r^}. There is no upper bound for tq. Let us take a fixed non-null value for tq. Then, 
can decrease continuously from = Tq to a lower bound given by 

(r^)min = rl- (S2)max , (34) 

where (s2)max can be calculated using a diagonal matrix K = diag(x+, x_, 0, ■ ■ ■ ,0). We have 
to maximize S2 = subjected to si = x+ + X- = const and > 0. We easily find 

that x+ = X- = Si/2 maximizes S2 which yields, using Eq. ( 126|) . 



~ 2{N-1 

In fact, we can write the two non-null eigenvalues x± of K as 



"^^^"^ - a% ^ J^-^ . (35) 



^± = V ^(/-l) ro±\/^-alrl. (36) 

Therefore, the SU(N) orbits are entirely specified by the value of r?^ = f^/rp G [otv,!]- In 
particular, = aj^VQ, instead of Eq. fl33|) . defines the maximally charge breaking space. 

The result f l35|) implies that, in the r^-space, the orbit space is restricted to a conical region 
between two coaxial cones: the lightcone, and the inner one defined by (l35l) . This observation 
might lead to non-trivial topological properties of the model, similar to what was described 
in [13]. 

Finally, when passing from the full to the neutral orbit space, one might be surprised that 
a single algebraic condition (!33|) reduces the dimensionality by 2N — 3 units, from 4{N — 1) 
(charge-breaking) to 2N — 1 (neutral). To show how this happens, let us write S2 in terms of 
doublets. We find 

S2{K)= Zab, where z^fe = (0i0a)(0l0fe) - (0l0b)(0l0a) • (37) 

l<a<6<Af 

There are A^(A^ — l)/2 quantities Zab, and each of them is non-negative due to the Schwarz 
lemma. Not all of Zab are independent, though. Suppose that all norms {<p\(pa) are fixed. 
Then, for the first three doublets, the quantities Zi2, z\z and Z2z are algebraically independent. 
However for any extra doublet, e.g. 0g, only two of z's, e.g. Zig, Z2q, can be chosen indepen- 
dently. Any further Zaq with a > 2 is not independent anymore but is linked to previous 2;'s 
by an algebraic relation, see a proof in Appendix [XI This is a consequence of the fact that 
we deal with doublets, not higher representations of the gauge group. Thus, for N doublets 
we have 2N — 3 independent Zab- Now, requiring that S2 = automatically sets all Zab = 0, 
which means that it is equivalent to 2A^ — 3 independent equalities. 



3.3 Root space 

Reparametrization transformation of the doublets, 4>a ^ 4>a = Uab(pb, described by a unitary 
matrix Uab G SU{N) corresponds in the adjoint space to a certain rotation of the vector r*: 
^ = OijTj, where Oij = Oij{U). The transformation matrix Oij belongs to the group 
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adj5't/(iV) (adjoint representation), which is only a proper subgroup of SO{N'^ — 1). It means 
that not all rotations in 5*0 (A^^ — 1) can be induced by reparametrization transformations. 

This fact restricts the way we can manipulate the adjoint orbit space. However, we always 
have a reparametrization freedom to bring any initial ii"-matrix to the diagonal form. In the 
adjoint space, it corresponds to certain allowed rotations of the entire orbit space that bring 
any point down to the (A^— l)-dimensional root space, which describes the diagonal i^-matrices. 
In the n space the N neutral orbits are represented by vertices of a regular (A^ — l)-simplex, 
while the charge-breaking orbit space is represented by the edges of this simplex, i.e. by the 
line segments joining the vertices. This gives the full description of the orbit space in the root 
space. For example, the orbit space restricted to the root space corresponds, for = 3, to the 
vertices and edges of an equilateral triangle while, for A^ = 4, it corresponds to the vertices 
and edges of a regular tetrahedron. The case A^ = 3 will be treated in more detail in Sec. HI 

The vectors in the root space can be parametrized in a very symmetric way in terms of A^ 
barycentric coordinates pi constrained by J2iLiPi = 1 Pi ^ 0- 



where en are the canonical matrices defined by {eij)ki = SikSji. Additionally, we can have at 
most two non-null pi, since rank K < 2. The vertices of the simplex, corresponding to the 
neutral orbit, are given by p = (1, 0, . . . , 0), (0, 1, . . . , 0), . . . , (0, . . . , 0, 1). Notice the matrices 
Qi are not all independent but obey Qi = 0. Various geometric features can be calculated 
explicitly by using the coordinates pi. 

The orbit space in the root space has a residual discrete symmetry with group Sn, related 
to the permutations of the doublets and corresponding permutations of the vertices of the 
simplex. Thanks to this freedom, we conclude that any neutral orbit can be brought to a 
predefined vertex, which means that all the points in the neutral orbit space are conjugate to 
each other, that is, can be mapped to each other by a reparametrization transformation. As 
for the charge-breaking points, one can always use the reparametrization freedom to place it 
on any predefined edge of the simplex, and even more, on any of the two symmetric halves 
of the edge. Therefore, if we are restricted to the points not conjugated by Sn, we get a 
line segment going from one vertex to the middle point of an edge. Such minimal set is in 
one-to-one correspondence to the SU{N) orbits in the fi-space. As we already discussed, one 
parameter can be chosen to characterize each point in the line segment, i.e., |n| G [a^, 1]. Then 
the whole orbit space V<i, can be recovered by SU{N) conjugation on this line segment and by 
varying tq. One can also recover the result fl5^ just from the shape of the orbit space in the 
root space by calculating the distance from the midpoint of an edge of the (A^ — l)-simplex to 
its center. 

3.4 Isotropy groups and SU{N) stratification 

Let us also describe the isotropy groups of the charge-breaking and neutral vacua, that is, the 
subgroups of the total reparametrization group that leave invariant a given point. 




(38) 



where we defined the traceless matrices 




(39) 
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Let us take a point in the neutral orbit space and bring it down to the root space, turning 
its iT-matrix into diag(f ^, 0, . . . , 0). It remains invariant under any U{N — 1) transformation 
that does not involve the first doublet as well as under a U{1) phase rotation of the first 
doublet alone. Since the bilinear are insensitive to the overall phase rotation, we get the 
isotropy group of the neutral vacuum SU{N — 1) x U{1), which is a (A^ — 1) ^-dimensional Lie 
group. 

Since the entire reparametrization group SU{N) has N"^ — 1 parameters, there are 2N — 2 
generators that do shift a chosen point along mutually orthogonal directions in the orbit space. 
Therefore, the neutral orbit space has 2N — 2 dimensions in the n space, and 2N — 1 dimensions 
in the r'^ space. This coincides with the calculations of Section 12.11 

Now take a generic point in the charge-breaking orbit space, with its i^T-matrix being 
diag(vf , 0, . . . , 0). It remains invariant under U{N — 2) transformation of the last N — 2 
doublets as well as phase rotations of the first two doublets. The isotropy group is therefore 
SU{N-2) X U{1) X f/(l), whose dimension is (A^- 2)^ + 1. The coset space SU{N)/{SU{N- 
2) X f/(l) X f/(l)) has dimension 4A^ — 6, which gives the dimensionality of the charge-breaking 
SU {N)- orbit, where the chosen point lies. Since we have a one-parametric family of such orbits 
by varying vf and w| but keeping the sum (r°) constant, we conclude that the dimension of 
the charge-breaking n-orbit space is AN — 5. In the r'^-space, the dimension is 4iV — 4, which 
again coincides with the counting of Section 12.11 

Now let us take a closer look at a point lying in the "maximally charge-breaking" or- 
bit space, that is, the one with i^'-matrix conjugate to diag(f^, f ^, 0, ... ,0). It corresponds 
to maximally charge breaking (MCB) orbit because S2, which quantifies charge breaking, is 
maximum for a fixed tq. It also corresponds to vectors n with smallest |fip possible and, 
therefore, lying on the surface of the inner cone. Such a point has a larger isotropy group than 
a generic charge-breaking point. Indeed, its isotropy group is now SU{N — 2) x SU{2) x U{1) 
of dimension (A^ — 2)^ + 3, therefore the dimension of the coset space (and of the maximally 
charge-breaking n-orbit space) is 4A^ — 8. When considering r^, it corresponds to a manifold 
of dimension 4N — 7. 

To summarize, we can group the SU{N)-oThits into classes of orbits according to its 
isotropy groups. A set of orbits with the same isotropy group is called a stratum [34] . For our 
problem, we have in general three strata for a fixed tq: 

(I) = 1, one (neutral) orbit, isotropy group SU{N — 1) U{1). 

(II) n'^ G {a%,l), continuous set of (charge-breaking) orbits, isotropy group SU{N — 2) 
f/(l)®f/(l). 

(III) fp = a^, one (maximally charge-breaking) orbit, isotropy group SU{N — 2) (g) SU{2) ® 
U{1). 

Notice that for A^ = 3 the strata I and III have the same isotropy group. 

As a final remark, we note that SU{N) / SU{N — 1), in fact, defines the space of A^-complex- 
dimensional vectors of unit absolute value (i.e. sphere S"^^"^). Its coset space with respect to 
the group U{1) of the overall phase rotations, {SU{N)/SU{N -1))/U{1) = SU{N)/{SU{N- 
1) X f/(l)), is by definition the complex projective space CP^~^. Thus, we recover the shape 
of the neutral orbit space just from its isotropy group. 
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4 Orbit space of 3HDM 



4.1 The three sets of coordinates 

In this Section we analyze the orbit space of the three-Higgs-doublet model in more detail. 

The Higgs field space of the 3HDM has 12 dimensions, hence the dimensionality of the 
charge-breaking and neutral orbit spaces are 8 and 5, respectively. They are embedded in the 
9-dimensional space of (ro,rj), i = 1, ... ,8, and are limited to the interior and the surface 
of the forward cone Tq — = 0. In the 8D space of "normalized" vectors n = r/vQ, the 
charge-breaking and neutral orbit spaces are 7D and 4D, respectively. 

The i^-matrix is a hermitean 3x3 matrix, which is decomposed via the unit matrix and 
the Gell-Mann matrices Aj, i = 1, . . . , 8: 



K = ro- i=l3 + r,A, . (40) 
The explicit expressions for the coordinates are 

ro = -^TtK = i= + 4<P, + 4<P,^ . (41) 

1^ rr^M ('^I'^i) - (<^2</'2) i4<Pi) + (0202) - 2(0503) 



2^^^"-'^' 2 ' 2^3 

ri = Re(0l02) , r2 = Im(0|02) , = Re(0l03) , 

rs = Im(0l03) , re = Re(403) , rj = Im(403) . (42) 

It is also useful to group the last six real coordinates (which we will refer to as the "transverse" 
coordinates) into three "complex coordinates": 

ri2 = ri + ir2 , r^^ = ri- ir^ , tq-j = rQ + ir-r . (43) 

The same indices accompany the normalized vectors n. 

The root space of the 3HDM is represented by the (71,3, n8)-plane (all the other rii = 0), 
shown in Fig. [H left. The neutral manifold intersects this plane by three distinct points P, 
P', P": 

P : is: oc diag(0, 0, 1) , ^3 = 0, rig = -1 ; 

P' : irocdiag(l,0,0), ^3 = ^, rig = ^ ; (44) 

/3 1 

P": ocdiag(0,l,0), 7x3 = rig = - . 

The charge-breaking manifold is represented by the three line segments joining these three 
points. Thus, the full orbit space in the root plane is given by the equilateral triangle (the 
2-simplex). Note that this triangle lies in the annular region between the circles of radii 1/2 
and 1, in compliance with (135|) . 

The triangle has the S3 symmetry, however the choice of coordinate used to describe it, 
and ns, breaks it. To restore this symmetry in the description, we introduce two additional 
coordinate sets on the same plane: (n3,?T,g) and (rig, rig), which are shown in Fig. [H right, by 
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Figure 1: (Left) The (ns, n8)-plane; all other = 0. Shown are the unit circle (section of 
the lightcone), the inner circle (dashed line), the three points P, P', P", from the neutral 
manifold, and the three line segments from the charge-breaking manifold. (Right) The same 
plane, but with three sets of coordinates: (ns, ng) shown in black solid lines, (^3, n'g) shown in 
dashed blue lines, and (n'g, n'g) shown in dash-dotted red lines. 



blue dashed and red dash-dotted axes. These coordinate sets are obtained from (723,77,8) by 
27r/3 and 47r/3 rotations, respectively: 

% = -2^3 ± , n^, % = T-^'^s - • (45) 

Each of the points P, P', P" can be associated with its "natural" coordinate set: 

P : ris = 0, ns = -I , P' : n'3 = 0, n'g = -1 , P" : n'^ = 0, n'^ = -1 . (46) 
The coordinates Ug, n'^, n'^ are closely related to the three barycentric coordinates 

l-2ng , l-2n' „ l-2n'l , „ ^ 

P = ^ , P = ^ , P = ^ , p + p + p = 1 , (47) 

which are proportional to the distances from a given point on the root plane to the three edges 
of the triangle. The three edges of the triangle, which describe the charge-breaking orbit space 
on the root plane, can be naturally parametrized by p = 0, p' = 0, and p" = 0. 

Thus, the points on the root plane can be described in a symmetric fashion using either 
{773,^3,^3} "with ns + ng + n'^ = 0, or {ns,ng,n'^} with ng + n'g + rz'g' = 0, or {p,p',p"} with 
p + p' + p" = 1. 

This symmetric description can be extended to the entire orbit space V<i>. Indeed, the 
27rA;/3-rotations on the root plane are generated by a cychc permutation of doublets: 

{0;, 0'2, 0'3} = {02, 03, 0l} , {01, 02, ^s) = {03, 01, M ■ (48) 

This permutation changes the "transverse" coordinates (143!) according to 

Wi2^ n'i5^ ^Gj} = {^67, ni2, ^45} , {n'/a, ^'5, n'^j} = {^45, ngy, ^12} . (49) 
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In other words, structures in the entire orbit space can be described in an exphcitly S^- 
symmetric way using coordinates 

p, p', p", ni2, ^12 5 where p + p+p" = l. (50) 

The unit sphere is represented in terms of these coordinates as 



2 



|2||/|2||//|2 1 /ciX 

+ F12I + F12I =1- (51) 



The i^-matrix of 3HDM also takes a very symmetric form: 



Kab = 4>l<Pa = ro 



ni2 V3p" n'* 



12 



(52) 



\ n'l* n'la VSp J 
Finally, the three quantities Zab (1371) can be written as 

2^12 = 3pV' - |ni2|^ > , zrs = "ipp - 1^12^ > 5 -^23 = 3p"p - |%2p > . (53) 

4.2 d-condition 

In 3IIDM, the if-matrix is a positive-semidefinite matrix with zero determinant, [10]; thus, 
the list of constraints on the coordinates of r'^ truncatesH at ( 128|) . which we will refer to as 
the "d-condition" . In the n-space this condition can be written as 

r , 3^2-1 
VMijkniUjnk = • (54) 

In order to select out the neutral manifold, we accompany the (i-condition with fp = 1, which 
makes it 

VsdijkniUjnk = 1 . (55) 

Alternatively, the neutral manifold can be defined even more compactly with the aid of the 
"star-product" (m *n)fe = y/SdykmiUj (V-product in Ref. [TT]): 

n^ = l, n * n = n . (56) 

Let us now write the (i-condition (15^ explicitly using the well-known values of dijk'- 

y/S / 2 : 2 2 2\ 3io ( 2 , 2, 2 -\- -\- n-j 



-6 • —[niniriQ + nin^rij - 77-2714 77.7 + 712715770) = • (57) 



^One can check explicitly that the higher order equations become identities. For example, thanks to the 
relationdijcdcki+djkcdca+dkicdcji = (SijSki+SjkSu+SkiSji) which holds for N = 3, we get r|^i/7'i?'j7-j,r; = r*, 
which makes S4{K) = satisfied automatically. 
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It can be rewritten in terms of symmetric coordinates (1501) : 

2 

p|ni2|^ + p'|%2|^ + p"|%2p - ^PP'p" - y|Re(ni2ra'i2<2) = , (58) 

which exposes the symmetry of the orbit space. One can also arrive at this expression 
directly from deiK = using representation (1521) for the K-matrix. 

For the neutral manifold, the ci-condition can be simplified further. Let us recall that the 
lightcone condition fl33|) implies that all three Zi2, Z23, 231 are equal to zero. Using fl53|) . and 
denoting the sum of the phases of ni2, ^12 7, one can cast the d-condition for the neutral 
orbit space into 

PpV'(1 - C0S7) = 0. (59) 



4.3 The local properties of the orbit space 

So far, we have described the shape of the orbit space on the root plane (713, ng), with all the 
transverse coordinates ni2 = n^^ = uqi = 0. Let us now gain an intuitive picture of how the 
orbit space extends into the transverse space. 




Figure 2: Generating a 4D cone from the triangle by mixing the first two doublets. 



In principle, the entire orbit space can be reconstructed by applying the full group of 
adjS'f/ (3) of orthogonal transformations of n induced by unitary SU{3) transformations among 
the doublets to the triangle on the root plane. To make this result more visual, let us first 
consider the subgroup of adjS'[/(3) induced by SU{2) transformations between the first two 
doublets: 



/ cosfe^T 
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The corresponding transformation of vectors rii 
point with coordinates 



(60) 





1/ 

fii brings a point on the root plane to the 



ni = — sin a cos(/3 — 7) 77,3 , 
ns = cos aris , ns = ns , 



n2 = — sin a sin(/3 — 7) 77,3 
n45 = = . 



(61) 



The SU{2) subgroup of such transformations, which we call 7^-rotations, applied to the triangle 
sends it to the surface of a 4D cone lying in the 71,45 = ^^67 = subspace with the apex at 
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point P, which is schematically illustrated by Fig. [2J Indeed, the upper edge of the triangle is 
mapped to the 3D ball 

3 1 

nj + nl + nl< ^s = ^45 = ne? = , (62) 

which serves as the base of the cone and which is nothing else but the orbit space of the 
2HDM. The two other edges of the triangle are mapped to the lateral surface of the cone 

nl + nl + nl=^^±pl. (63) 

Note that the neutral orbit space is represented in this 4D cone by the apex and by the "rim" 
of the based, the sphere + rig + ^3 = 1 at = 1/2. 

The similar constructions arise from mixing of other pairs of the doublets. Namely, the 
SU{2) subgroup that mixes 02 and 03 (T^'-rotations) keeps point P' invariant and sends the 
triangle to a 4D cone with the base 

< + ri? + <<^, n', = ^, (64) 

which lies in the subspace ni2 = ^^45 = 0. Finally, mixing (pi and 03 (7?."-rotations) generates 
a similar cone with apex at P" lying in the subspace ni2 = = 0. 

Thus, in very loose terms, the shape of the orbit space can be described as follows: it is 
a manifold stretched between differently oriented three 4D cones. However, when visualizing 
this picture, one should remember that in fact there is no distinction between the base, the 
lateral surfaces of these cones and of the part of the orbit space that is stretched between the 
cones. The neutral orbit space, being nothing but CP , is a homogeneous manifold, so it looks 
the same at every point. For the charge-breaking orbit space, rather similarly, there is a "flat 
face" going through each point. Some additional hints for visualization of CP^ are given in 
[351 [28]. 

To make these observations more precise, let us calculate the sectional curvatures along all 
mutually orthogonal directions at any point in the orbit space. 

We start with a point located on a charge-breaking manifold. By an appropriate reparametriza- 
tion transformation we bring it to the root plane and place it, for example, on the upper side 
of the triangle, where its position is described by = 1/2 and some n^. We know that this 
point lies inside a fiat 3D ball. Hence, there are three directions (parallel to axes rii, n2 and 
n^), along which the orbit space is fiat in the vicinity of the selected point. 

We are left with four other directions, along n^, n^, ng, and ny. One can shift into these 
directions by performing small rotations introduced above. Note that from the point of view 
of the 71'- and 7^"-rotations, the upper edge of the triangle is located at the lateral edge of 
the corresponding cone, which brings in some curvature. 

Explicitly, let us apply to a point on the upper edge the sequence of an 7?.'-rotation with 
an infinitesimal a' and an 7^"-rotation with an infinitesimal a", all the other angles 7', 
7" being arbitrary (the order of the two transformations is irrelevant for this calculation). We 
get shifts of 71,45 and ngy linear in the small angles 



18 



and shifts in 77,3, ng which are quadratic in small angles: 



^ 16 
16 



2 \ / 9 



9 \ (9 



(66) 



For the charge-breaking manifold, where is a flat direction, we need to keep track only of 
the changes in ng, for which we get: 

~ - , (67) 

where the curvature radii along directions uq, uj (-Re, 7) and along directions 724, (-^4,5), are 

Thus, the charge-breaking orbit space has locally the shape of an ellipsoidal cylinder, with 
three fiat directions and two pairs of curved directions with sectional curvature radii i?4^5 and 

-^6,7- 

We now repeat this calculation for a point at the neutral manifold, for example, point P. 
We again perform two infinitesimal rotations TZ' and TZ" and calculate shifts of the coordinates. 
This time we must take care of shifts of all eight coordinates nj. These rotations give linear 
shifts in small a' and a" to the four transverse coordinates, 

a/3 ^/S ^/S ^/S 

and quadratic shifts to the other coordinates 

Sn,, ^ ^a'a" e^^^""^') , Sn^ ^ ^{a"^ - «'') , 5ns ^ ^{a'^ + a"^) . (70) 
The overall quadratic shift is 

5n = ^{5n,Y + {dn^f + {5n^y + {dngf ^ ^{a'^ + a"^) ^ I^^^^P + I'^ngrP _ ^^^^ 

4 V 3 

Thus, the curvature radius of the neutral manifold is Rq = a/3/2 regardless of the direction of 
the shift. Since this holds true at every point of the neutral manifold, it means that the neutral 
manifold is an example of spherical space forms (a manifold of constant sectional curvature). 
This comes as no surprise: it is known that an even-dimensional spherical space must be a 
sphere or a complex projective space |36] . 

Note that the curvature radius Rq does not and should not coincide with the largest 
curvature radius of the charge-breaking manifold near the rim. In loose terms, the fact that 
the neutral orbit space is located on the unit sphere gives to the neutral points more curvature 
with respect to the adjacent charge-breaking points. 
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4.4 Duality property of the orbit space 

The orbit space of 3HDM has an additional duality property, which does not hold for a generic 
A^: if a ray along direction n goes through the neutral orbit space, then a ray in the opposite 
direction, — n, points towards a maximally charge-breaking point. Since the charge-breaking 
points lie on the sphere |n| = 1/2, we find that the maximally charge-breaking orbit space is 
homothetic to the neutral orbit space with the scale factor of 1/2. 

This property can be easily proved in the root plane: if the neutral point is characterized 
by the i^-matrix diag(0,0,f^) = y(l — y/SXs) {\n\ = 1), then the opposite point corresponds 

to the i^-matrix diag(f ^, f ^, 0) = ^(1 -|- ^As) (|n| = 1/2). This property clearly depends on 
the number of the diagonal elements and does not generalize for higher A^. However, at iV = 4 
another observation can be made: if n points towards a maximally charge-breaking point, 
then so does —n. For example K = diag(l, 1, 0, 0) is opposite to = diag(0, 0, 1, 1) with the 
same tq. That is, the maximally charge-breaking orbit space in the four-Higgs-doublet model 
is centrally symmetric. 

5 Conclusion 

In this paper we initiated an analysis of the general A^-Higgs-doublet model. Focusing only 
on the scalar sector of the model, we considered here a specific question: how to efficiently 
describe the space of gauge- invariant bilinears of Higgs fields in NHDM (the orbit space). We 
characterized the orbit space as a certain algebraic manifold embedded in the Euclidean space 
and studied some of its algebraic and geometric properties. The general construction was 
illustrated with the case of = 3, for which more detailed calculations were presented. 

For general NHDMs for N > 2, compared to the N = 2 case, there arises a general 
and distinct feature of the orbit space: the orbit space is no longer convex, i.e., for two 
arbitrary points x'^ {Ki) and y'^ {K2) in V$ [M^(A^;2)], the line segment joining them may 
not be entirely contained in V$ [M^(iV;2)] [21]. For example, for Ki = diag(f^, f ^, 0) and 
K2 = diag(0, u^, V?), their middle point is \{Ki + K2) = |diag(f ^, f ^ + u^, v?) which no longer 
has rank 2 or smaller. The exception is the case of two neutral points, as explained in Sec. 13. II 
More particularly, we showed there is a "hole" in the orbit space of constant tq, such that in 
n-space it is constrained inside the annular region of radius \n\ = 1 (lightcone) and |n| = a^v 
(inner cone). In other words, for Tq > 0, we can not reach |r| < aj^VQ. This feature will bring 
very distinct possibilities to the symmetry breaking patterns of the potential as well as to the 
conditions for bounded below potentials. Some of its consequences will be further detailed in 
a forthcoming work [25] . 

We also commented on a remarkable similarity between the orbit space of NHDM and the 
state space of an A^-qudit in quantum information theory. We sketched a small "dictionary" 
between some objects in these two branches of theoretical physics, and we think that this link 
should be explored further. 

The next step of this analysis, the study of the NHDM Higgs potential and its symmetries, 
is done in the companion paper [25] . That study is also conducted in the orbit space and uses 
many of the results of the present paper. We hope that the methods presented in these papers 
will boost systematic exploration of the wealth of structures hidden in the general NHDM. 

It is clear that a very similar mathematics arises not only in multi-doublet models, but 
also in models with A^ copies of Higgs fields in other representations (scalars, triplets, etc). It 
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is therefore conceivable that even more comphcated Higgs sectors can be treated along these 
lines. Other possible applications could be found in the condensed matter physics, where 
group- invariant potentials depending on several interacting order parameters are often used, 
[37] . An example where the methods of 2HDM were used to understand the general Ginzburg- 
Landau model with two order parameters can be found in [38]. 
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A Relations among Zat 

In order to show that among A^(A^ — l)/2 quantities Zab in the NHDM there are only 2A^ — 3 
algebraically independent quantities, we need to prove the following statement. Take any 
four doublets, e.g. 0i to 04 with known norms, suppose that Zi2, Zi^, z^, Z23, are 

also known (we assume here a generic situation when all Zat are non-zero). Then 2:34 is not 
independent and can take at most two different values. If all of Zu, zi^, Z14, Z23, -224 happen to 
be zeros, then Z34 = 0. 

We first note that the doublets are vectors in the space C^. Therefore, if Z12 7^ 0, all the 
doublets can be decomposed in the basis of 0i and 02 : 

03 = Cl0l + C202 , 04 = C?101 + C?202 • (72) 

Note that the absolute value of the scalar product (0l02) is known, |(0i02)P = (0j0i)(0202) — 
Zi2, its phase 6*12 is not. Let us now introduce the "vector" product of two doublets: 

[0„ X 0,] ^ 0,+0O - 000+ , (73) 

where superscripts + and refer to the upper and lower components of the doublets. One can 
check that 

[0a X 0fe]*[0, X 0b] = (0l0a)(0|;0b) - (0l0fe)(0l0a) = Zab • (74) 

This leads to 

I , |2 ^24 I , |2 ^14 /„p.N 

— , = — , \d2\ = — . (75) 

Zl2 Zi2 Z\2 

Therefore, decomposition (172]) turns into 

— e^''^0i + ./^e^''^02, 04 = ./^e^«Vi + ./^e^«^02. (76) 

2^12 V ^12 V ^12 V ^12 



1 2 _ ^23 






|C2| 


^ Z\2 ' 
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The phase differences rj2 — rji and ^2 ~ are both related to the (unknown) phase 612'- 



2:12 (0404) 



2:23(0101) + ^13(0202) + 2^^132:23 I (0102) I COs(6'i2 + ^2 - Vl) , 



^24(010l) + 2;i4(0202) + 2^-214^241(0102)1 COs(6'i2 + 6 " 6) • 



(77) 



Now, the quantity 2:34 can be written as 



^12^34 = ^23^14 + ^24^13 - 2^/^23^14^24^13 COs(?7i + 6 " ^2 - ^l) • 



(78) 



But 



COs(r7i + 6 - ^2 - ^l) = COs[(6'i2 + 6 - 6) - (^12 + V2- Vi)] 

cos(6'i2 + 6 - 6) cos(6'i2 + ?72 - ?7i) ± I sin(6'i2 + 6 - 6) sin(6'i2 + m - Vi) \ , 

which can be expressed in terms of known cosines of 6*12 + ^2—^1 and ^12 + ri2 — t]i- This 
proves an algebraic relation between 2:34 and the other quantities without the need to know 
6^12- The sign ambiguity here means that two different values of 2:34 can result. However, if 
-2^23-2^14-2:24-2:13 = 0, then 2:34 is uniquely determined. 

Note that if all of Zi2, 2:13, 2:14, 2:23, ^24 happen to be zero, it means that all four doublets are 
proportional to each other, and therefore, 234 must be zero as well. 

As a remark, let us analyze in more generality the phenomenon of multidimensional re- 
duction imposed by a single condition fl5^ . A more general situation can be envisaged. The 
space of N X N hermitean matrices with rank equal or lower than r < N, which we can 
denote by Mh{N;r), has dimension r{2N — r). To define Mh{N;r), we need Sn{K) = 0, 
r < n < N. Despite only one constraint Sr-i{K) = being further required to restrict 
Mh{N;r) to Mh{N;r — 1), the dimensionality is indeed reduced by 2(A^ — r) + 1. In our 
case, we have r = 2 and the amount of dimensional reduction from Mh{N; 2) to Mh{N; 1) is 
exactly 2A^ — 3. Therefore, any single condition Sr-i{K) = necessary to restrict Mh{N;r) 
to Mh{N;r — 1) should contain multiple independent conditions in the same way S2{K) = 
is equivalent to various conditions Zab = 0, as proved in this appendix. 

B Maximal set of gauge invariants 

We will show here how we can choose a maximal set of algebraically independent gauge 
invariants 0^0a = Kab, corresponding to the iN — 4 degrees of freedom of the doublets 0a. 
If all bilinears 0^0a, a,b = 1, . . . , N, a < b, were functionally independent, A^^ real parameters 
would be necessary for parametrization. 

Firstly, we should use the fact that a general non-null K matrix ([6]) has rank two or one. If 
it has rank two, it is always possible to choose a set of two linearly independent lines (columns) 
of as a basis of the space spanned by all the A^ lines (columns); otherwise only one line 
is linearly independent and this case can be treated easily. By appropriately labeling the 
doublets we can choose the first and second lines to be non-null and non-parallel. In that 
case, since A' is a hermitean matrix, we can choose the set in Eq. iQ as the minimal set of 
gauge invariants. It is easy to see that they can be parametrized by AN — 4 real parameters, 
considering that Ku, K22 are real. We should assume Ku 7^ and A'22 7^ because, e.g., the 
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case Kii = directly implies Kia = Kai = 0, a > 1. Such property follows directly from the 
fields language but it also can be thought as a consequence of 

5^|i^la|'< ^11 ( 5^ i^aa), (79) 

which follows from the Schwarz inequality. Thus any null diagonal element implies an entire 
null line and column of K. 

It remains to be shown that all other Kab-, with a, 6 > 3, can be written entirely in terms of 
the set in Eq. ([9]). Let us show how to calculate the elements in the third line. The calculation 
of any other element follows analogously. By hypothesis, we can write any element in the 
third line as a linear combination of the corresponding element in the first and second lines: 

fTsa = aKia + PK2a , a > 3 . (80) 
But the the same coefficients a,P relate the elements in the first and second columns as 

Ks2 = K;^ = aK^2 + fiK22. (81) 
Equations flHTl) can be rewritten as 

(7^31 i^32) = (« /3) Kf^ , (82) 

where K^^^ is a 2 x 2 submatrix (minor) of K containing only the elements Kab, with a = i,j 



^3) 



and b = Thus we can invert equation (!82|) to obtain 

= [K,, K,2) {Kt^V (J;;) , 
or 

det(iriJ)ir3a = (i^3i K,2) adj(irg)) {^^^^ , (84) 

where adj denotes the adjoint matrix. Notice det(ii'];2 ) = -2^12 is non-null by hypothesis. 

We can rewrite Eq. (1831) in a more compact form if we define the 2-dimensional complex 
vector 

xl={K,a K2a) , a = l,...,iV. (85) 
Then any matrix element Kab can be calculated as 

Kba = A{Kf^)-\a. (86) 

Surprisingly, the expression in Eq. (1861) is valid not only for a, 6 > 3, but for all a, 6 = 1, . . . , iV. 
However, for a,b = 1,2, it leads to trivial identities. 

Equations (!83|) and (186!) are direct consequences from the fact that any 3x3 submatrix 
of K has null determinant for rankK < 2. For example, Eq. (18^ is equivalent to calculate 

the determinant of a 3 x 3 matrix constructed with the blocks KI2 ,Xa,X3^ ^3a by cofactor 
expansion along the third column. 
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One can identify the role of the coefficients a, (3 in Eq. (180 p if we recognize the equation as 
the expansion 

03 = a0i + /302 , (87) 
contracted to (p^. Hence, the coefficients of the hnear expansion 

(pa = Cal4>l + Ca2/302 , (88) 

are solutions of 

{Kal Ka2) = (Cal C.s) K^^ ■ (89) 

The hermitean conjugate of Eq. ( 189|) can be also written 

Xa = C*iXl + C*2X2 • (90) 

If ranki^r=l, we would have Kba = KbiKia/ Ku, Ku ^ 0, for all a, 6 = 1, . . . , A^. 



C Characterization of SU{N) orbits 

The vector space spanned by the N x N hermitean matrices, containing K, is isomorphic to 
R^^, where the vectors r'^ live. The mapping between these spaces were given by Eq. ( fTOl) 
and it is valid even if we generalize i^' to be a general N x N hermitean matrix. The action 
of the reparametrization group SU{N) on K is defined by Eq. ([7]). Such action divides the 
space into SU{N) orbits. Each of these orbits can be uniquely characterized by a set of 

SU{N) invariants functions Sk{K), k = 1, . . . , N, defined in Eq. fl2T]) . Therefore, any orbit 
can be represented by a point in one connected region of a A^- dimensional diagram whose 
axes represent Sk- There is only one connected region because we can vary the eigenvalues 
continuously, keeping, for instance, a decreasing order. 

The reparametrization group action in Eq. ([7]), however, defines naturally two invariant 
spaces (irreducible representations) which allows the splitting 

K = Ko + K, (91) 

where Kq = si{K)l/N and K = K — Kq = TiXi is the traceless part of K. Hence, K is 
the component of K that transforms non-trivially under SU{N) while Kq is an invariant. 
This means that the invariants Sk{K), k > 2, are not fundamental but have contributions of 
the trivial part Kq, already in si{K). We can use Sk{K), instead of Sk{K), for k > 2, which 
obviously are invariant and does not depend on Kq. Let us denote Sk = Sk{K) and Sk = Sk{K). 
The relation between the two sets {sk} and {sk}, k = 2, . . . , N, can be obtained by comparing 
Eq. ([19]) to 

N 

det(Al -K) = det(Al - K) = + ^(-l)^SfcA^-^ , (92) 

k=2 

where X = X — ^. The relation between Sp and Sp, for p > 3, is 

^p-sp - -I:r^'')(^)V.+(P-l)©(^)^ (93) 

fc=i 

- i:m(|)V. +©(!)'. (94) 

k=l 
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At last, all invariants Sk = Sk{K) can be calculated using Eq. f l2ip and written in terms of 



|lr[(r,\r] = r(2,...,„r,r,,---r. 



(95) 



where the tensors F^") were defined in Eq. (!30|) . For example, 



§2 = -2Tr[(rjAi)^] = -f^ , 
S3 = |Tr[(riAi)^] = dijknrjrk , 



(96) 
(97) 
(98) 



All Sk can be written in terms of the terms of Eq. (195 p with equal or lower order. 

It is important to notice that for general N x N hermitean matrices K, not restricted 
to positive semidefinite rank two matrices, the characterization of the SU{N) orbits would 
involve more than one invariant, besides tq. For instance, for a value of f^, there would be 
infinitely many distinct orbits that have to be further characterized by higher order invariants. 
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